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ABSTRACT 

We identified 24 Si IV absorption systems with z < 1 from a blind survey of 49 low-redshift 
quasars with archival Hubble Space Telescope ultraviolet spectra. We relied solely on the char- 
acteristic wavelength separation of the doublet to automatically detect candidates. After visual 
inspection, we defined a sample of 20 definite (group G = 1) and 4 "highly-likely" (G = 2) 
doublets with rest equivalent widths W r for both lines detected at > 2>aw T - The absorber line 
density of the G = 1 doublets was cWsirv/dA = 1.4±g;| for logiV(Si +3 ) > 12.9. The best-fit 
power law to the G = 1 frequency distribution of column densities /(iV(Si +3 )) had normalization 
k = (1.2±g;|) x lCT 14 cm 2 and slope a N = — 1.6±g;|. Using the power-law model of /(V(Si+ 3 )), 
we measured the Si +3 mass density relative to the critical density: flgi+ 3 = (3-Tti 7) x 10~ 8 for 
13 < log iV(Si +3 ) < 15. From Monte Carlo sampling of the distributions, we estimated our value 
to be a factor of 4.8^t 3 'g higher than the 2 < z < 4.5 (f2gj+3). From a simple linear fit to n Si +3 
over the age of the Universe, we estimated a slow and steady increase from z = 5.5 — > with 
dft si+ 3 /dt agc = (0.61 ± 0.23) x 10~ 8 Gyr -1 . We compared our ionic ratios A(Si +3 )/A r (C+ 3 ) to a 
2 < z < 4.5 sample and concluded, from survival analysis, that the two populations are similar, 
with median (iV(Si +3 )/iV(C +3 )} = 0.16. 
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1. Introduction 

The signatures of the cosmic enrichment cy- 
cle are etched into the processed gas and re- 
flected in its metallicity, elemental abundances, 
density, and/or spatial distribution. Measuring 
these quantities con strain models of ga lactic feed- 
back processes fsee lBertschingerlll998l and refer- 
ences therein). Currently, quasar absorption-line 
(QAL) spectroscopy is the best tool for probing 
the IGM. 
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The ultraviolet (UV) transition Si IV AA1393.76, 
1402.77 A is a well-studied doublet in QAL sur- 
veys covering 1.5 < z < 5.5. From the observ- 
ability perspective, the Si IV absorption lines are 
valuable for the following three reasons. First, 
Si IV absorbers can be observed outside the Lycv 
forest since they have rest wavelengths red-ward 
of Lycv. A1215. This reduces the effect of blending. 
Second, they are observable from ground-based 
telescopes when they redshift into the optical 
passband at z > 1.5. Third, they constitute a 
doublet with characteristic rest wavelength sep- 
aration (9 A) and equivalent width ratio (2:1, 
respectively, in the unsaturated regime). When 
these distinctive criteria are met, we can be fairly 
confident that the pair of absorption lines are a 
Si IV doublet. 

From the astrophysics perspective, the Si IV 
doublet is a strong, observable transition of silicon. 
The abundance of silicon is predominately driven 
by Type II supernovae at z > 1, with an increas- 



1 



ing fraction from feedback from asymptotic giant 
branch stars more recently (jOppenheimer fc Dave 



20081 : 

icon traces 



C IV systems with column den s ities of 10 



Wiersma et all 120091 ) . Thus, at z > 1, sil- 
most abundant metal. 



to 10 14 cm~ 2 (Kim et al 



2003 : Aguirre et al.l 12004 ). Indeed, any variation 



2002; iBoksenberg et al 



oxygen, the 

Oxygen itself is difficult to study since its strong 
transitions are blue- ward of Lya (e.g., O VI 
AA1031, 1037), if not also blue- ward of the Lyman 
limit (e.g., O IV A787 ) . By using Si I V absorption 
as a tracer of oxygen, ISongailal (|2001l ) constrained 
the IGM metallicity to be > 10~ 3 5 Z Q at z = 5. 
Therefore, the fraction of cosmic star formation 
that occurred before z = 5, or within 1 Gyr of the 
Big Bang, was > 10~ 35 . 

ISongailal ( 2001 ) also measured Sl Si +3 , the Si +3 
mass density relative to the critical density, for 
2 < z < 5.5. The Si +3 mass density may represent 
the largest contributio n to the silicon mass density 
fisi at these redshifts (lSongailall200ll) but still con 



in the ionic ratio iV(Si +3 )/iV(C +3 ) may be dom- 
inated by the variation in the me tallicity of the 
absorbing gas ((Bolton fc ~vlell2010h . However, ev- 
idence for both a break and strong evolution in the 
UVB have been detected in some studies of Si IV 



and C IV absorbers (|Songailal ll998. 2005) 



The gas giving rise to 12 < logiV(Si +3 ) < 14 
absorption has T < 10 4 9 K and n early constant 
rSi/C l » +0.77 for 1.5 < z < 4.5 ((Aguirre et all 
2004)0 In simulations of the IGM, most (possibly 
all) of the silicon is located in distinct clouds of 
metal-enriched gas. At z ~ 3, the clouds have 



radii « 1 Mpc (proper; IScannapieco et al" 2006 ) 



stitute only a small fraction of fisi ((Aguirre et al 
20041 ). For z = 4.5 -> 2, O si +3 was roughly 
constant at « 1.2 x 10~ 8 for absorbers with col- 
umn densities 13 < logN(Si +3 ) < 150 For 
z = 5.5 — > 4.5, the Si +3 mass density may have 
increased by an order of magnitude. Subsequent 
studies have largely supp orted these broad trends 
in r?s;+3 (non ) evolution 



and could be considered filamentary structure. 
However, by z — 0, the enriched clouds are ac- 
tually th e extended gaseous halos fa 100 kpc) of 
galaxies ((Dave fc Oppenheimerll2007l ). Absorbers 
in the circum-galactic medium would likely be 
subjected to a softer ionizing background (due 
to the increased stellar contribution) than the 
Haardt fc Madaul (|l996l ) UVB typically used in 



(Bokscnberg et al.l 2003; 



Songailall2005l ; IScannapieco et al.ll2006l ). 

Observations of gas bearing both Si IV and 
C IV AA1548,1550A absorption offer constraints 
on the shape, spatial extent, and/or evolution 
of the ionizing ultraviolet background (UVB; 



IGM studies. Local sources (i.e., star-forming 
galaxies with a non-zero escape fraction of ionizing 
photons) may be the most important contributor 
to the background. If the background for Si IV 
absor bers were softer at z < 1, [Si/C] would be 
lower ((Aguirre et al.ll2004j ). 

The current work finalizes our analysis of 



Boksenberg et al . 2003; I Aguirre et al.l2004l IScannapieco etaafcl iival Hubble Space Telescope (HST) UV spec 



2006T ). The ionization threshold for Si 



++ 



-to-Si+ 3 

is 2.5 Ryd and for C++-to-C+ 3 , 3.5 Ryd. Thus, 
the ionic ratio N (Si +3 ) / N (C +3 ) is affected by the 
shape of the UVB at these energies, and the fluc- 
tuations in the ratio spatially and in time could 
constrain the patchiness and evolution, respec- 
tively, of the UVB. 

Of particular interest is the effect of He II reion- 
ization on the UVB at z 3, which is expected to 
boos t Si IV absorption and s uppress C IV absorp- 
tion ( Madau fc Haardt 2009). Several studies find 
no evidence for a sharp break in the shape of the 
UVB at z w 3 or even significant evolution in its 
shape for z « 4.5 —> 1.5 from studies of Si IV and 



tra, gathered prior to Servicing Mission 4 (UT 
July 2009). This is the largest survey for Si IV 
systems at z < 1 to date and covers the last eight- 
billion years of the cosmic enrichment cycle. In 
other words, the net effect of cosmic star formation 
(and feedback) at the 'end' (i.e., z = 0) is con- 
strained by observations of Si IV absorption at low 
redshift. Also, the survey of z < 1 Si IV absorbers 
provides a baseline for similar, high-redshift sur- 
veys, which are currently more numerous. The 
data reduction and analysis metho ds used in this 



paper are described in detail in iCooksev et al 



( 2010L hereafter Paper I), to which the interested 
reader is referred. 



5 We adjust quantities from other studies to our adopted 
cosmology: Ho = 70kms _1 Mpc -1 , £7m = 0.3, and Qa = 
0.7. 



6 We adopt the following notation: [Si/C] = log(ns;/rtc) — 
log(nsi,Q/nc,©), where nx is the volume density of ele- 
ment X. 
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This paper is organized as follows: we present 
the data processing and sample selection in Sec- 
tions [2] and El respectively; we analyze and dis- 
cuss the frequency distribution, its moments, and 
V(Si +3 )/A(C+ 3 ) in Section |U and Section [5] is a 
summary. 

2. Data, Reduction, and Measurements 

We conducted a blind survey for z < 1 Si IV 
systems in the Hubble Space Telescope (HST) UV 
spectra of 49 low-redshift quasars, which makes 
the current work the largest low-redshift Si IV 
study to date. We included spectra from the 
Space Telescope Imaging Spectrograph (STIS; 
pre-Servicing Mission 4) and the Goddard High- 
Resolution Spectrograph (GHRS). The STIS 
cchclle spectra, taken with the E230M grating, 
provided most of the search path length (see Fig- 
ure!]}, but we also searched the other STIS echelle 
grating (E140M) and the GHRS echelle (ECH-B) 
and long-slit (G160M, G200M, G270M) gratings. 
Spectra from the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) covered the transitions with rest 
wavelengths X r < 1100 A (e.g., higher-order H I 
Lyman lines) . All spectra had resolution with full- 
width at half-maximum FWHM < lSkms" 1 and 
signal-to- noise ratio S/N > 2pix _1 . 

The spectra were retrieved from the Multi- 
mission Archive at Space Telescope (MAST)0 
The reduction and co-addition of multiple ob- 
servations foll o wed t he algorithms described in 
Cooksev et al. (2008). The spectra were normal- 
ized semi-automatically. All reduced, co-added, 
and normalized spectra are available online, even 
those not explicitly searched in this paper H 

The central wavelength (and redshift) of an ab- 
sorption line was measured by the optical depth- 
weighted mean of the pixels within the wavelength 
bounds (Xi,Xh) defining the absorption line (see 
Table [l}. The rest equivalent widths W r were 
measured with simple boxcar summation, and the 
column densities (e.g., A(Si +3 )) were measured 
with the apparent optical depth method (AODM; 



7 See http://archive.stsci.edu/ 

8 See 



http: / /www. ck.org/~xavier /HSTSilV/ 
normalized spectra, the continuum fits, the Si 1 



for the 
V candi- 



date lists, and the Monte Carlo completeness limits for all 
sightlines as well as the completeness test results for the 
full data sample. 



ISavage fc Sembachl 1991 ). The doublet column 
density A(Si +3 ) is either: the variance-weighted 
mean of the measurements of both lines; the col- 
umn density from the one line with a measure- 
ment; the greater lower limit; or the mean, when 
the limits of the lines define a finite range. 

The co-moving path length AX sensitive to 
Si IV doublets with W r > io~w r in both lines 
was estimated from Monte Carlo simulations. For 
these simulations, we replaced all automatically- 
detected features in the archival spectra with ran- 
dom noise drawn from surrounding pixels. For 
each redshift bin (8z = 0.005) of each spec- 
trum, we distributed 10 3 Si IV doublets with a 
range of column densities, Doppler parameters, 
and number of components. Then we measured 
the iV(Si +3 ) and W r . 1393 limit per bin at which 
95% of the doublets were automatically recov- 
ered (see Figure [2]) . The available path length in- 
creases with increasing A r (Si +3 ) and W r ,i393; for 
the strongest absorbers (logiV(Si +3 ) > 14.1 and 
^,1393 > 174mA), Al = 18. 

3. Si IV Sample Selection 

Absorption lines with observed equivalent 
width W bs > 3aw obs were automatically detected 
in the spectra. Candidate Si IV doublets were 
identified based solely on the characteristic wave- 
length separation of the doublet (w 1940 km s -1 ). 
First, each Wobs > ^ a w oba line was assumed to be 
Si IV 1393, and any automatically-detected line 
that was near the location of would-be Si IV 1402 
was adopted as such. If no Wobs > 3o"iv obs line 
existed, an upper limit was set on Wobs from the 
spectrum. Second, any automatically-detected 
line not already tagged as Si IV 1393 or 1402 was 
assumed to be the latter, and an upper l imit was 
set on Wobs,i393 from the spectrum (see iPaper 1 . 
for details). A summary of all candidate Si IV 
doublets are given in Table Q] 

Other common absorption lines (e.g., Lya and 
C IV) were associated with the candidate sys- 
tems in a similar manner. For example, an 
automatically-detected absorption line would be 
identified as a candidate Lya line if it had ob- 
served wavelength A ~ A Q (1 + z C and), where A Q is 
the rest wavelength of Lya and z can( j is the red- 
shift of the candidate Si IV doublet. The FUSE 
spectra were useful in this step, since they cov- 
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Table 1 

Si IV CANDIDATES SUMMARY 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


^1393 


5t>abs 


X T 


A; 


Xh 






log AT 


Clog N 


Flag 




(kms -1 ) 


(A) 


(A) 


(A) 


(mA) 


(mA) 
















MRK335 (z cm = 


0.026) 










-0.00001 




1393.76 


1393.40 


1394.04 


212 


15 


13.55 


0.04 


495 




-6.9 


1402.77 


1402.41 


1403.06 


129 


15 


13.54 


0.05 




0.00643 




1393.76 


1402.41 


1403.03 


125 


14 


13.23 


0.05 


317 




-7.3 


1402.77 


1411.48 


1412.10 


< 32 




< 13.15 a 






0.01078 




1393.76 


1408.67 


1408.91 


< 20 




< 12.65 a 




173 




0.7 


1402.77 


1417.78 


1418.02 


37 


12 


12.97 


0.14 










PG0117+213 (zom 


= 1.493) 










0.63473 




1393.76 


2278.01 


2278.98 


622 


70 


> 13.84" 




258 




10.0 


1402.77 


2292.74 


2293.72 


< 83 




< 13.59 






0.63639 




1393.76 


2280.25 


2281.14 


< 116 




< 13.41" 




160 




-0.5 


1402.77 


2295.00 


2295.90 


240 


41 


> 13.72 a 






0.63962 




1393.76 


2284.12 


2286.44 


478 


67 


> 13.95 




259 




94.0 


1402.77 


2298.89 


2301.23 


< 167 




< 13.81 












TONS210 (z om = 


0.116) 










-0.00070 




1393.76 


1392.56 


1393.09 


81 


12 


13.04 


0.06 


367 




-4.4 


1402.77 


1401.57 


1402.10 


< 25 




< 12.77 






-0.00002 




1393.76 


1393.34 


1394.13 


272 


14 


13.63 


0.03 


431 




7.8 


1402.77 


1402.35 


1403.15 


110 


17 


13.47 


0.06 




0.00129 




1393.76 


1395.39 


1395.69 


29 


9 


12.55 


0.14 


301 




-8.0 


1402.77 


1404.42 


1404.72 


< 26 




< 12.79 







a log N measured by assuming W r results from the linear portion of the COG. 

Note. — Summary of Si IV doublet candidates by target and rcdshift of Si IV 1393. Upper limits arc 2a 
limits for both W r and log AT. The binary flag is described in Scction[3l 

(This table is available in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 
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Fig. 1. — Schematic of redshift coverage for the 
current survey. The number of spectra with cov- 
erage of the Si IV doublet is shown as a function 
of Z1393 (histogram). The STIS E140M spectra 
covered 21393 < 0.21. The STIS E230M spectra 
typically covered 0.6 < 21393 < 1.2. The redshift 
range 0.21 < 2:1393 < 0.6 was covered by some 
E230M spectra as well as STIS medium-resolution 
gratings and GHRS spectra. The redshifts of the 
doublets detected with W r > 3<Jw T hi both lines 
are shown with the hashes across the top. The top 
and middle rows indicate the redshifts of the 18 
unsaturated and two saturated doublets, respec- 
tively, in the definite group (G = 1). The bot- 
tom row shows the redshift of the four unsaturated 
"highly-likely" (G = 2) doublets. 

ered the Ly^ A1025, O VI, and C III A977 lines 
for z < 0.1 candidate systems. The presence of 
associated lines increased the confidence of our 
identifications. 

Each candidate doublet was assigned a machine- 
generated, binary flag that scaled with the number 
of desired characteristics of a true Si IV doublet. 
The characteristics (flags) are as follows: 

256: W r ,l393 > 3<7h/.,1393; 

128: W r .uo2 > 3(7^1402; 

64: the ratio ^^393 : is in the range 

of 1 : 1 to 2 : 1 plus/minus the propagated 
error of the ratio; 

32: the optical depth- weighted centroids of the 
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o r ./.,....,....,.. 1 

50 100 150 

Wr,1393 

Fig. 2.— Redshift path length AX(log 7V(Si +3 )) 
and AX(W rj i393) as a function of Si IV 1393 col- 
umn density (top) and rest equivalent width re- 
spectively (bottom). These estimates are based 
on Monte Carlo analysis and correspond to 95% 
completeness limits. 

Si IV lines have <5v a bs| < 10kms _1 H 
16: there exists a candidate Lya line with 

W r ,a > 3(JW r ,a] 

8: the 1393 line is outside of the Lya forest 

4: and outside the H 2 forest; 

2: the smoothed AOD per pixel of the doublet 
lines agree within la for > 68.3% of the 
pixel; and 



The wavelength separation between a line at z a t, s and a 
Si IV line at 21393 is <5f a bs = e(z abs - 21393)/(1 + Z1393)- 
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1: there exists one or more candidate lines (not 
H I) associated with the candidate doublet 
and with W r > 3aw r ■ 

All candidates were visually inspected by at 
least one author, and the candidates with both 
lines detected at W r > 3aw r were reviewed by 
two or more. We agreed upon 22 definite Si IV 
systems, which constitute the "G = 1" group (see 
Figure [3]). Of these, 20 have both doublet lines 
detected with rest equivalent width W r > 3aw r 
and constitute the group on which we based our 
conclusions. The G = 1 sample has a median 
redshift (z) = 0.906, log7V(Si+ 3 ) > 12.9, and 
W r ,i393 > 66 mA. 

We also defined a small, "highly-likely" (G = 
2) sample of six systems (see Figure 0]). These 
doublets are typically found in regions with low 
S/N and/or do not have other lines associated with 
them, which would increase the confidence of our 
identification. The four of these with both lines 
detected at > 3aw r were combined with the G 
= 1 sample for some analyses and identified in 
tables and figures as G = 1+2. Details of all (G = 
1+2) absorption systems are given in Table[2] The 
properties of the all Si IV doublets are summarized 
in Table [3] 

All doublets are more than 1000 km s 1 out- 
side of the Galaxy and 3000 km s -1 blue- ward of 
the background quasar. There were no Si IV ab- 
sorbers without associated Lya absorption, when 
the spectral coverage existed to detect Lya. We 
combined Si IV doublets into one absorption sys- 
tem when their optical depth-weighted centroids 
had 5w a bs < 250 km s -1 . 

The observed absorber line density is the sum 
of the number of absorbers, each weighted by the 
path length sensitive to their iV(Si +3 ) or Wi-,1393- 
For the G = 1 sample, <W SiW /dX = l-4±g:| 
(dA/~siiy /dz = 2. 7±glg) for logA^(Si +3 ) > 12.9. 

As in IPaperl . we conducted Monte Carlo sim- 
ulations to measure the rate that pairs of z < 1.5 
H I Lyman forest lines satisfy the characteristics 
of Si IV doublets and were potentially included in 
our sample as doublets. The contamination rate 
of forest lines masquerading as Si IV was small, 
less than 5% of dAfsnv/dX or an expected lj^ 
false doublet. If any doublet in our sample were 
false, it would be one without other associated ab- 
sorption lines. There are four such doublets, all in 



the G = 2 sample. Our expectation that forest 
lines could mimic Si IV doublets drove us to de- 
fine the "highly-likely" G = 2 group. Though we 
provide the results for analyses of the G = 1+2 
sample, we only discuss the results from the G = 
1 sample and base our conclusions on that; so we 
concern ourselves no further with the effects of the 
H I Lyman forest contamination. 

3.1. Comparison with Previous Studies 

There have been three recent surveys for Si IV 
systems at z < 1 using at least some of t h e HS T 
spectra analyzed here: Milutinovic et al.l (|2007l ): 
iDanforth fe Shulll (|2008l ). and IPaper J . Here we 
briefly compare our blind doublet search results 
with these other studies. In the first two, they 
identified Lya lines first and then sought associ- 
ated transitions 

identified 17 Si 



Milutinovic et al 



(12007ft identified 17 Si IV 
doublets in the eight STIS E230M spectra that 
they surveyed. We independently identified 13 
of those absorbers. The remaining four doublets 
were explicitly listed as questionable by these au- 
thors. We identified the G = 2, Z1393 = 0.65536 
absorber in the PG01 17+213 spectrum, one o f 
the eight surveyed by iMilutinovic et al. ( 2007 ). 
though they did not detect it. Since only the dou- 
blet was detected, they would not have found it 
with their Lya-targeted search. 



Danforth fc Shulll (|2008l ) surveyed all of the 



STIS E140M spectra for many transitions, includ- 
ing Si IV. They did not require that both doublet 
lines be detected with W r > 3aw r ■ They found 20 
Si IV doublets for a line density dJVsnv/dX = bt\ 
for W r > 30 ml. 

In contrast, we found seven Si IV absorbers in 
the E140M spectra. There were two (G = 1) dou- 
blets th at we identified in our S i IV-targeted sur- 
vey but lDanforth fc Shulll (|2008l ) did not find: the 
z 1393 = 0.00 5 72 doublet in QSO-123050+011522 
and the z 1393 = 0.13846 one in PG1116+215 . 
As mentioned in lPaper A [Danforth fc Shulll (|2008h 
missed the former doublet because the Ly/3 line in 
the FUSE spectra was suspect. The latter doublet 
has W rt i402 < 3(7^,1402, which might explain why 
they did not identify it. 

Of their 20 systems, we agree with five and 
independently identified them in our Si IV- 
targeted survey. The remaining 15 doublets from 
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Danforth fc Shulll ([20081 ) were not included in our 
sample for at least one of the following reasons: 
one or both lines detected were at W r < 3aw r (12 
doublets); one line was blended with a Galactic 
line (1); it was a system intrinsic to the back- 
ground quasar (1); and/or the doublet was ob- 
served in the Oth order of E140M (2), which was 
exclude in our reduction. 

We recovered a ll Si IV doublets that we previ- 
ously identified in [Paper II because of the associa- 
tion with a C IV system. We also determined that 
the G = 1, zi548 = 0.240 1 C IV doublet from 
Paper J is actually O I A1302 and Si II A1304 as- 
sociated with the G = 1, Z1548 = 0.47436 system, 
which corresponds to our G = 1 Si IV system with 

^1393 = ^1548- 

4. Analysis 

4.1. Frequency Distributions 

Analogous to the luminosity function used in 
galaxy studies, observers define the column den- 
sity frequency distribution /(iV(Si +3 )) to be the 
number A7V of Si IV doublets per column density 
interval AiV(Si +3 ) per path length AX (see Figure 
[5]). With a maximum likelihood fitting algorithm, 
we fit a power law to /(iV(Si +3 )) as follows: 



PG01 17+213 (01:20:17.25) 



Lycc Forest 



z„„= 1.04801 



/(X(Si+ 3 )) 



/ JV(Si+ 3 ) 
V N 



(1) 



where k is the normalization with unit cm 2 ; Nq = 
10 13 ' 5 cm -2 ; and a/v is the slope. We often refer 
to ku, which is the normalization k scaled up by a 
factor of 10 14 . The frequency distribution was fit 
over the range k>g-/V m i n = 12.84 to log iV max = 15, 
with special treatment of the saturated absorbers 
(log iV sat = 14) in the maximum likelihood anal- 
ysis (see Paper J ). Briefly, the number of dou- 
blets with 14 < log7V(Si +3 ) < 15 (i.e., saturated) 
was a constraint in our likelihood function. The 
best-fit parameters were: k u — l.lSj^ggCm 2 and 



a N - -l.bl_ 031 . 

The fit parameters for /(iV(Si+ 3 )) and /(W^.,1393) 
(discussed below) are given in Table SJ We esti- 
mated the 68.3% confidence limits (c.l.) in the 
power-law normalization and slope by tracing a 
contour where SC = InC — ln£ max = —1.15 on 
the likelihood surface C, which included 68.3% of 
its area. Then, the 68.3% c.l. (what we will loosely 
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Fig. 3.— Velocity plots of 22 G = 1 Si IV sys- 
tems. The regions of spectra around each absorp- 
tion line are aligned in velocity space with respect 
to the rest wavelength of the transition and Z1393. 
Saturated transitions are indicated with the (red) 
'S'; any transition with column density detected 
at < 3 cr are noted with a (red) 'W.' Systems with 
z i393 < A Q (1 + z cm ) are labeled with a (red) "Lya 
Forest." The regions used to measure W r and 
logiV(Si +3 ) are shown by the dark outline. The 
flux at zero and unity are shown with the dash-dot 
lines (blue and green, respectively); the vertical 
dashed line indicates v = Okms -1 , corresponding 
to the optical depth-weighted velocity centroid of 
the Si IV 1393 transition. (The G = 1 velocity 
plots are available in their entirety in the online 
journal. One figure is shown here for guidance re- 
garding the form and content.) 



refer to as "1-<t errors," hereafter) were defined as 
the difference between the k and a extrema on the 
contour and the £ max values. Errors in quantities 
derived from the frequency distributions (e.g., the 
Si +3 mass density discussed below) are estimated 
in a similar fashion. The likelihood surface is not 
Gaussian, and the "2- and 3-er errors" are de- 
fined by SC = -3.15 (95.4% c.l.) and -6.1 (99.7% 
c.l.), respectively. Therefore, the larger confidence 
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Table 2 
Si IV SYSTEMS SUMMARY 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


21393 


<5l>abs 


A r 








(mA) 


logiV 


Clog N 


G 


Flag 




(kms -1 ) 


(A) 


(A) 


(A) 


(mA) 











PG0117+213 (Zem = 1.493) 



0.65536 




1393.76 


2306.70 


2307.89 


273 


28 


> 13.73 




2 


482 




9.2 


1402.77 


2321.62 


2322.82 


122 


24 


13.54 


0.09 






1.04801 




1393.76 


2853.75 


2854.89 


194 


18 


13.54 


0.07 


1 


471 




-10.5 


1402.77 


2872.21 


2873.36 


83 


14 


13.37 


0.07 








-111.3 


1215.67 


2486.96 


2490.84 


1530 


21 


> 14.88 










-27.2 


1206.50 


2469.68 


2471.49 


446 


19 


> 13.67 










-26.2 


1260.42 


2580.19 


2581.61 


145 


11 


13.16 


0.03 












B0312-770 (2 




0.223) 










0.20253 




1393.76 


1674.90 


1676.82 


652 


78 


> 14.16 




1 


511 




7.6 


1402.77 


1685.73 


1687.67 


543 


57 


> 14.40 










-15.0 


1215.67 


1460.66 


1463.09 


1677 


20 


> 14.92 










-7.3 


1025.72 


1232.52 


1234.29 


1139 


20 


> 15.59 










11.2 


977.02 


1174.06 


1175.78 


912 


67 


> 14.44 










51.6 


1031.93 


1240.64 


1241.62 


496 


17 


> 14.94 










58.2 


1037.62 


1247.48 


1248.47 


310 


14 


14.87 


0.02 








21.7 


1206.50 


1449.82 


1451.60 


763 


19 


> 13.94 










61.5 


1238.82 


1489.57 


1490.45 


137 


7 


13.86 


0.02 








62.1 


1242.80 


1494.36 


1495.24 


94 


9 


13.99 


0.04 








39.1 


1260.42 


1515.44 


1516.45 


490 


12 


> 13.92 














PKS0405-123 (z om = 


= 0.573) 










0.16712 




1393.76 


1626.40 


1626.95 


117 


15 


13.24 


0.06 


1 


499 




5.4 


1402.77 


1636.92 


1637.47 


103 


23 


13.49 


0.10 








-18.1 


1215.67 


1417.80 


1419.38 


876 


19 


> 14.60 










-0.9 


1025.72 


1196.77 


1197.50 


444 


20 


> 15.19 










-22.2 


977.02 


1139.53 


1140.73 


500 


22 


> 14.26 










12.1 


989.80 


1154.85 


1155.66 


214 


12 


> 14.64 










-30.1 


1031.93 


1203.74 


1204.67 


391 


24 


> 14.77 










-50.7 


1037.62 


1210.37 


1211.31 


218 


37 


14.80 


0.11 








2.3 


1206.50 


1407.81 


1408.40 


236 


11 


> 13.38 










-6.4 


1238.82 


1445.53 


1446.12 


108 


15 


13.81 


0.06 








-7.4 


1242.80 


1450.17 


1450.76 


76 


13 


13.94 


0.07 








-3.3 


1260.42 


1470.73 


1471.28 


148 


12 


> 13.28 









Note. — Si IV systems by target and rcdshift of Si IV 1393. Upper limits are la limits for both W r and logiV. 
The column densities were measured by the AODM, unless "COG" is indicated, in which case, the limit is from 
assuming W r results from the linear portion of the COG. The definite Si IV doublets arc labeled group G — 1, 
while the "highly-likely" doublets arc G — 2. The binary flag is described in Scction[3l 

(This table is available in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 
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Table 3 
Si IV DOUBLET SUMMARY 



(1) 

Target 


(2) 
G 


(3) 

21393 


(4) 

WV,1393 

(mA) 


(5) 

W r ,1402 

(mA) 




(6) 

log JV1393 




(7) 

log JVi4 2 




(8) 

log AT(Si+ 3 ) 


PG0117+213 


2 


0.65536 


273 ± 28 


122 ± 24 


> 


13 


.73 




13. 


.54 ±0 


.09 


13.54 ±0.09 




1 


1.04801 


194 ± 18 


83 ±14 




13. 


.54±0.07 




13. 


37±0 


.07 


13.43 ±0.05 


B0312-770 


1 


0.20253 


652 ± 78 


543 ± 57 


> 


14 


.16 


> 


14 


.40 




> 14.40 


PKS0405-123 


1 


0.16712 


117 ± 15 


103 ±23 




13 


.24 ±0.06 




13 


.49 ±0 


.10 


13.27±0.05 


PKS0454-22 


1 


0.47436 


635 ± 52 


475 ± 56 


> 


13 


,84 a 


> 


14 


.02" 




> 14.02 




1 


0.48331 


112 ±37 


< 61 


> 13 


Of) 11 


< 


13. 


.34 




[13.09, 13.38] 


HE0515-4414 


2 


0.68625 


249 ± 18 


68 ± 17 


< 


13 


,71 b 




13 


.36 ±0 


.10 


13.36±0.10 




1 


0.94047 


233 ± 13 


56 ±11 


< 


13 


.58 b 




13. 


.15 ±0 


.08 


13.15±0.08 




1 


1.14760 


565 ± 22 


343 ± 29 




13 


.98 ±0.03 




13 


.98 ±0 


.04 


13.98 ±0.02 


HS0810+2554 


2 


0.55455 


250 ± 54 


< 160 


> 


13 


,44 a 


< 


14 


.06 




[13.44, 13.76] 




2 


0.83143 


262 ± 30 


271 ±31 




13 


.67 ±0.08 




13. 


91±0 


.06 


13.77 ±0.05 


MARK132 


2 


0.93166 


73 ± 12 


86 ± 12 




12 


.95 ±0.07 


< 


13 


31 b 




12.95 ±0.07 


PG1116+215 


1 


0.13846 


34 ±9 


< 18 




12 


.68±0.10 


< 


12 


.64 




12.68 ±0.10 


PG1206+459 


1 


0.92690 


777 ± 19 


402 ± 29 




14 


.08±0.01 




14 


.06 ±0 


.03 


14.08 ±0.01 




1 


0.93429 


143 ± 11 


98 ± 11 


> 


13 


.51 




13. 


.60 ±0 


.06 


13.60 ±0.06 


PG1211+143 


1 


0.05118 


66 ±4 


44 ±5 




12 


.96 ±0.03 




13. 


.04 ±0 


.05 


12.97 ±0.03 


MR.K205 


1 


0.00428 


122 ± 11 


75 ± 13 




13 


.25 ± 0.04 




13. 


.33 ±0 


.07 


13.27 ±0.04 


QSO-123050+011522 


1 


0.00572 


66 ±8 


63 ±9 




12 


.96 ±0.05 




13 


.21 ±0 


.06 


13.02 ±0.04 


PG1248+401 


1 


0.77305 


472 ± 19 


348 ± 17 


> 


14 


.01 




14 


.06 ±0 


.02 


14.06 ±0.02 




1 


0.85485 


211 ± 21 


138 ±31 




13. 


.49 ±0.05 




13. 


.59 ±0 


.10 


13.50 ±0.04 


PG1630+377 


1 


0.91432 


123 ± 15 


98 ± 12 




13 


.26 ±0.06 




13 


.47 ±0 


.06 


13.33 ±0.04 




1 


0.95279 


202 ± 15 


136 ±21 




13 


.49 ±0.04 




13 


.59 ±0 


.07 


13.51 ±0.03 


PG1634+706 


1 


0.65351 


145 ± 11 


71 ±9 




13. 


.30 ±0.04 




13. 


.25 ±0 


.06 


13.29 ±0.03 




2 


0.81813 


< 15 


21 ±7 


< 


12 


.25 




12 


.71 ±0 


.13 


12.71 ±0.13 




1 


0.90560 


182 ± 4 


32 ±5 


< 


13 


,68 b 




12 


.92 ±0 


.07 


12.92 ±0.07 




1 


0.99035 


240 ±7 


162 ±6 




13 


.73 ±0.02 




13 


.74 ±0 


.02 


13.74 ±0.01 




1 


1.04106 


167 ± 5 


96 ±8 




13 


.32 ±0.02 




13 


.35 ±0 


.03 


13.33 ±0.01 


PHL1811 


1 


0.08094 


108 ±7 


74 ±9 


> 13 


.44 




13 


.45 ±0 


.05 


13.45 ±0.05 



'MogiV measured by assuming W r results from the linear portion of the COG. 
b Limit due to blended line. 



Note. — Summary of Si IV doublets by target and redshift of Si IV 1393. The definite Si IV doublets are labeled group G 
— 1, while the "highly-likely" doublets are G — 2. Upper limits arc la limits for both W r and log TV. The adopted column 
density for the Si IV doublets arc listed in the last column (sec Scction[2t. 
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PG01 17+213 (01:20:17.25) 



Lycc Forest 



z„„ = 0.65536 
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0.5 



PL, 
— 




i 1.2 



o 
Z 



SilV 1393 s 




0.5 



-100 100 

Relative Velocity (km s~ l ) 

Fig. 4. — Velocity plots of six G = 2 Si IV sys- 
tems. (See Figure for description of velocity 
plot.) (The G = 2 velocity plots are available in 
their entirety in the online journal. One figure is 
shown here for guidance regarding the form and 
content.) 



limits on the best-fit power-law parameters are, 
formally, as follows: 2a kli = +0.79/ — 0.54 cm 2 ; 
2a aN = +0.45/ -0.53; 3a kli = +1.18/ - 0.70 cm 2 ; 
and 3cr QN = +0.62/ - 0.76. 

There was no observed break in /(V(Si+ 3 )), 
and no break has been observed at high redshift. 
There must be a break in order to limit the number 
and mass of Si IV absorbers to finite quantities. 

We have measured a slope consistent with those 
from high -redshift studie s, though potentially 
shallower. ISongailal (|l997l ) measured a^ = —1.8 
z < 3 and -2 for 3 < z < 3.54 

(2006) 



for 2.16 < 



Songailal (|200ll ) and IScannapieco et al. 
stated that ajy = —1.8 matched their observed 
/(V(Si+ 3 )) well, which covered 1.78 < z < 5.29 
and 1.5 < z < 3.1, resp ectively. From thei r Lya- 



targeted z < 0.4 survey, iDanforth fc Shulll (|2008h 
measured a N = —1.92 ± 0.17 for a sample with 
different selection criteria and from a survey with 
different methodology. 



The definition of the equivalent width fre- 
quency distribution f (W r ,13m) is similar to that of 
/(V(Si +3 )), and it was also fit well with a power 
law, with W r fi = 150 mA, JV^min = 57 mA, and 
W> ;max = 796 mA. The W r limits reflect the ex- 
trema of the observed values, ±lcr. The best-fit 
values for /(W+,1393) were: k 3 = 3.481^2, mA 
and aw = —1.28 



+0.47 
-0.47- 



The W r k-i is the nor- 



malization scaled up by a factor of 10 3 . The 
larger confidence limits on the best-fit power- 
law parameters for f(W r ,i3Q3) are, formally, as 
follows: 2a k3 = +2.55/- 1.73 mA _1 ; 2a aw = 
+0.79/ -0.80; 3a ks = +3.81/ - 2.21 mA" 1 ; and 



3(T 



-1.10/-1.13. 



4.2. Si IV Absorber Line Density 

We measured the observed Si IV line density by 
the sum of the number of absorbers, each weighted 
by the path length sensitive to their V(Si +3 ) or 
W r 1393 (see Figure [2]). For the G = 1 sample, 
<W SiW /dX = lAtVa (cLVsiiv/dz = 2.4±g;J) for 
logV(Si+ 3 ) > 12.9. 

bmce /(V(Si +3 )) is modeled well by a power 
law, we can integrate Equation [T] to estimate the 
Si IV absorber line density dJVsnv /dX for a given 
column density limit Vi; m : 



IV 



(7V(Si+ 3 ) > V lim ) 



dX 



Tyl+ajv 
lim 



1 + CtN 



N* N 



(2) 

This is useful for comparing to high-redshift stud- 
ies, which typically do not match our observational 
limits. In Figure|6l we show aWs; rv / dX as a func- 
tion of Z1393 for Vii, n = 10 13 cm~ 2 . 

The Si IV line density does not increase sig- 
nificantly from z ks 3 — > 0. Our integrated, 
z < 1 'cbV Si Tv/dV = 1.2±°'j 
is co nsistent within la of the 
from IScannapieco et al. (120061) 



(logV lim = 13) 
(z) = 1.9 value 
and 2a of their 



(z) — 2.7 one. The formal errors on the inte- 
grated dA/sirv/dV are 2ad7 V '/dx = +1-3/— 0.6 and 
ia^M/dx — +1.2/ — 0.8, derived from the aN and 
fci4 errors discussed in Section [4.1l There is no ev- 
idence of evolution in dAfsnv / dX from z = 1 — > 
(i.e., within our sample). 

We estimated the high-redshif t line densities 
and er rors in the following manner. IScannapieco et al 



(|2006l ) published their /(V(Si +3 )), from which we 



estimated the power-law normalization fci4 given 
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Table 4 

Si IV FREQUENCY DISTRIBUTIONS SUMMARY 



(1) 

G 


(2) 

(z) 


(3) 
zi 


(4) 


(5) 
AA 


(6) 
Limits 


(7) 

<LVsiiv/dz 


(8) 

dA/"siiv/dX 


0) 

fi Si+3 X 10 8 


(10) 

fe a 


(11) 

a 


(12) 
Pks 




Column Density 


1 


0.90560 
0.91432 


0.00428 


1.14760 


20 
18 


(12.84, 15.00) 
(12.92, 14.40) 


2 7+ 1 - 1 

Z - '-0.8 

2 4+ ' 4 




i 9 + 0.5 
1 " i -0.4 
1 4+0-2 


r> *7-l +2.82 
°- 1 1 -1.68 

> 1.63 


1 1s +0.45 
1 - 18 -0.36 


-1 61+"- 2f! 


0.526 


1+2 


0.85485 
0.90560 


0.00428 


1.14760 


24 
22 


(12.84, 15.00) 
(12.92, 14.40) 


3 

°- J --0.9 

2 9 +0 - 4 




i r + 0.5 

■"■•°-o.4 
! 7 +o.2 

1 - ' -0.2 


, ,,+2.81 
*- u -l.T7 

> 1.94 


1 43+ - 48 


l.DO_ 2 8 


0.400 


Equivalent Width 


1 

1+2 


0.90560 
0.85485 


0.00428 
0.00428 


1.14760 
1.14760 


20 
24 


(57, 796) 
(66, 777) 
(57, 796) 
(66, 777) 


3 0+ 12 

O r, + 0.3 

z ' z -0.3 
Q r + 1.3 
°'°-1.0 

2 6 +0 3 

ZD -0.3 




1 4+0-5 

i O + 0.2 
i,o -0.1 

1.6+21 

— 0.5 
1 r + 0.2 
1 '°-0.2 




3.48ti;ti 
4.21+1;^ 




0.533 
0.563 



a The power-law coefficient & has units of 10 14 cm 2 for the column density section and 10 3 mA 1 for the equivalent width section. 

Note. — Parameters from the maximum likelihood analysis for f(x) — k (x / xq)™ , where x — iV(Si +3 ) (or W r ,i393) anc ^ is AT — 10 13 ' 5 cm - 2 
(W r ,o — 150 mA). For each G — 1 or 1+2 subsample, the first row summarizes the maximum likelihood analysis and the second row, the observed 
quantities. dAf Si i V /dX , listed in the first subsample row, is the integral of /(7V(Si +3 )) (/(W r( i393)) from log7V(Si +3 ) — 13 to infinity (W r ,i393 — 
50 mA to 796 mA) with the best-fit k and a. Also in the first subsample row, the integrated dJVsilv/d^ = °Wsi Yy/dX ■ dXfdz, where the latter 
term is evaluated at (z) . The observed dA/snv/d^ an d dA/siiv/dX are from the sum of the total number of doublets, weighted by the path length 
available to detect the doublet, based on its N(Si+ 3 ) or W r 1393- ^gj+3; listed in the first subsample row in the column density section, is the 
integral of /(7V(Si +3 )) ■ 7V(Si +3 ) from 13 < logiV(Si +3 ) < 15 with the best-fit k and a. The observed Si+3 were from the sum of the unsaturated 
doublets, as given by J\f. P K q is the significance of the one-sided Kolmogorov-Smirnov statistic of the best-fit power law. 
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-1.5 
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-3.5 

-4.0 
1 



13.0 13.5 14.0 14.5 15.0 

logMSr 3 ) 



G= 1: cc w = -1.28 „ 
k 3 = 3.48 mA- 



G = 1+2: a„ 



-1.30 



k"=4.2\ mA 



*3- 
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log ff ri i393 (m A) 



Fig. 5. — Column density (top) and rest equivalent 
width (bottom) frequency distributions for G = 1 
and 1+2 samples (black and gray, respectively). 
The best-fit power laws, from the maximum like- 
lihood analysis, are the dashed lines, with slope a 
and normalization fci4 as given. The column den- 
sity saturation limit is indicated by the vertical, 
dotted line (top). 



the best-fit slope ozn- We measured f(N(Si +3 )) 
at several 7V(Si +3 ) and used the scatter as an es- 
timate of the errors. We a ssumed a^v = —1.8, 
mapieco et al. ( 20061 ) did, and estimated 
1.3 ± 0.9 cm 2 for (z) = 1.9 and fe 14 = 
0.5 ± 0.1cm 2 for (z) = 2.7. 



as 
ku 



4.3. Si +3 Mass Density 

The Si +3 mass density h as been measured by 
sever al 1.5 <z< 5.5 studies (ISongailal ll997l 1200 li 
2005 ; Scannapieco et all 120061 ). Typically, they 



0.0 0.5 1.0 



1.5 2.0 

Z 1393 



2.5 3.0 3.5 



Fig. 6.— Redshift evolution of dJ\f Siw /dX. The 
line density of Si IV absorbers has not increased 
significantly since z w 3. The (orange) stars 
are our z 
10 13 cm" 2 



< 



+3\ 



> 



1 integrated values for iV(Si 
The high-redshift values are inte- 
grated dA/siiv/d^f, with the slope and normal- 
iz ation estimated from the f(N(Si +3 )) published 



Scannapieco et al. I (|2006l black squares) 



measured the mass density, relative to the critical 
density p r .,a, by summ ing the detected absorbers 
( Lanzetta et al. 199lh : 



H a m c \ - 

"si+ 3 - — — 2^ 



iV(Si+ 3 ) 



(3) 



where the Hubble constant is Hq = 70 km s _1 Mpc - 
the mass of a silicon atom is ms; = 5 x 10~ 23 g; c is 
the speed of light; and p Ci o = 9.26 x 10~ 30 gcm~ 3 . 
The summed ri Si +3 approximates the mass den- 
sity in all absorbers with column densities within 
the observed ./V(Si +3 ) range, which is 12 < 
logiV(Si +3 ) < 15 for the high-redshift studies. 

Since we measured column densities with the 
AODM, we only have lower limits on A^(Si +3 ) for 
the strong, saturated doublets, which dominate 
. Thus, summing the column densities only 
resulted in a lower limit Ogi+s > 1.55 x 10 -8 for the 
16 unsaturated G = 1 doublets with log N (Si +3 ) > 
13. We exclude the two saturated doublets since 
they might have logA^(Si +3 ) > 15, which is the 
upper limit we have chosen for comparing f2g;+3. 
Including the two saturated doublets, results in 
the summed fi Si +3 > 2.3 x 10~ 8 for the 18 G = 1 
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doublets with log7V(Si +3 ) > 13. 

In order to compare with the high-redshift 
studies, we assumed the power-law formalism for 
/(iV(Si +3 )) and integrated the column density- 
"weighted" /(7V(Si+ 3 )) (i.e., its first moment): 

gome k ( N%£»-N3[°> 
cp c ,o 2 + a N \ JVq n 

(4) 

From the best-fit values for the G = 1 / (iV(Si +3 )), 



sample: ** uc 



n; 



The median of this distribution 



9. 



Si+ 3 



Si+ 3 



(3.71 



+2.82\ 

-i.es) 



x 10~ 8 for logJV„ 



13 and 



logiV max = 15. Since there has been no observed 
break in /(iV(Si +3 )), the column density limits are 
crucial to defining a finite Si +3 and comparing 
between surveys. 

The formal errors on the integrated i7 Si +3 are 
2a n = +5.54/ -2.36 and 3ct = +8.99/ -2.81. 
derived from the a at and ku errors discussed in 
Section O 

We plot the evolution of r2 si +3 as a function of 
the age of the Universe i agc in Figu re E The me- 
dian^ of the 2 < z < 4.5 studies (|Songailall200ll 



IScannapieco et al" 20061 whose values have been 
adjuste d to mat ch our iV(Si +3 ) limits and cosmol- 



ogy, see lPaoer if l is (r! Si+ 3) 2 ^ 4 - 5 = (0.77i°; 3 i) x 
10~ 8 , which is shown by the (blue) lines in Figure 
(left). We have detected, with > 99.8% confi- 
dence, an increase in fisi+ 3 horn high-to-low red- 
shift. The z < 1 is a factor of 4.8l 3 g higher than 
the 2 < z < 4.5 median. 

The confidence limits on the median high- 
redshift value and the increase in fl Si +3 were es- 
timated based on Monte Carlo sampling of the 
distributions. First, we drew 10 6 realizations of 
the high-redshift data sets fi 2 ^- 4 ' 5 , assuming 
they had log-normal errors. Then, we measured 
the median of each set: (fij^f- 4 ' 5 )^. The median 
Si +3 mass density quoted above ((f2 si +3) 2 - z - 4 - 5 ) 
was the median of these (10 6 ) (Sljyj|?- 4 ' 5 )i values. 

To Monte Carlo sample our z < 1 measurement, 
we used the likelihood surface discussed in Section 
14.11 and the Metropolis-Hastings algorithm to ap- 
propriately sample the k-ajy parameter space 10 6 
times. For each of these random pairs, we com- 
puted the integrated Si +3 mass density (see Equa- 
tion , resulting in a low-redshift Monte Carlo 



10 When we quote median values with "l-<7 errors," these "er- 
rors" are actually the difference between the median (i.e., 
50th percentile) and the values at the 15.9th and 84.2nd 
percentiles. 



was: (fisi+a)^ 1 = (3.70t\H) x 10~ 8 E3 The ratio 
of the low- to high-redshift Monte Carlo samples 
(i.e., (^m~c~ 45 )'/^mci) * s a distribution where 
the median is 4.8^ 3 'g and the ratio is greater than 
unity at the 99.8% c.l. 

A least-squares minimization of a linear model 
to over t are for the 2 < z < 5.5 obs ervations 

(|Songailall20"0lT IScannapieco et ai~ll2006l) and our 
z<l value yielded: dfl Si + 3 /di age = (0.61+0.23) x 
10~ 8 Gyr, as shown by the (red) lines in Figure [7] 
(right). This toy-model slope agrees well (< 2cr) 
withd^c+s /dt ago = (0.42 ± 0.2) x 10~ 8 Gyr from 
Paper J for the equiv alent redshift sample bu t for 
C IV absorbers, and iD'Odorico et all (|2010h de- 
tected a smooth increase in flc+3 from z — 2.5 — > 
1.5, wh ich reaches the z < 1 9l c +3 measured in 



Paper J . A linear fit to only the 2 < z < 5.5 obser- 
vations indicated no temporal evolution in risi+ 3 
(i.e., dft Si+ 3 /dt ag c = (0.10 + 0.56) x 10- 8 Gyr). 

Modeling f2gi+ 3 as evolving linearly with time 
was not a physically motivated exercise but one 
method to evaluate whether our 9l si +3 indicated a 
significant increase compared to the high-redshift 
observations. Since the inclusion of our z < 1 
value resulted in a statistically significant rate of 
increase for f2gi+ 3 > we likely have detected a true 
increase of the Si +3 mass density at low redshift, 
though proof must await a larger z < 1 survey. 

Our results indicate that any increase in Si +3 
at z < 1 is likely due to an increase in the number 
of high-column density absorbers (i.e., shallower 
o>n compared to high redshift), since d/Vsiiv/d^ 
is nearly constant from z « 3 — > 0. In general, 
the Si +3 mass density is dominated by the high- 
column density absorbers, so even a small increase 
in their frequency will significantly change Slgi+ 3 ■ 

4.4. 7V(Si +3 )/iV(C+ 3 ) 

As mentioned previously, the ionic ratio iV(Si+ 3 )/ 
N(C +3 ) has been used to study the shape and/ 
or evolution of the UVB. In order to construct a 
complete sample of systems with coverage of both 
doublets, we measured the upper limit for N(C +3 ) 
(or iV(Si +3 )) when the doublet was not detected 



11 This value is in excellent agreement with our integrated 
value for the best-fit k and ajv and indicates a smaller 
spread, within the quoted confidence limits, than we for- 
mally adopted (see Table |4}. 
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Fig. 7. — Time evolution of Si mass density relative to the critical density. The int egrated for 
13 < log A r (Si +3 ) < 15 for the G = 1 sample, is the (orange) star. The value from iDanforth fc Shull 
([20081 ) is the (green) triangle an d is not an inde pendent measureme nt of We compare d our z < 1 

" (|2006l . black, open 



measurement only to those from ISongailal (|2001 
square). The median, and its "errors" (see Section |4~3 



black crosses) and iScannapieco et al. 

of the 2 < z < 4.5 values are shown on the left with 



the (blue) solid and dashed lines, respectively. The simple linear fit to risi+ 3 over t agc is shown on the right 
by the (red ) solid line, whe re the dashed lin es are the \-a ra nge of the fit. For reference, we also show the 
values from lSongaila ( 1997 . gray circles) and Songaila ( 2005 , gray, filled squares); the latter of which uses a 
pixel-optical depth method to measure fisi+ 3 - 



in as sociation with the targeted Si IV (or C IV, 
from iPaper J ) doublet but the spectral coverage 
existed. Due to how the S/N changed throughout 
any spectrum, a Si IV-targeted survey was not 
sufficient to define a complete C IV sample and 
visa versa. Ultimately, there were 12 detections 
and 12 lower limits for z < 1 A(Si +3 )/A(C+ 3 ), 
with log7V(Si+ 3 ) > 11.9 and logV(C+ 3 ) > 13.37. 

In Figure we com pare our 7V(Si +3 )/ 7 V(C+ 3 ) 
sample with that from Boksenberg et all ([2003). 
We reproduced their Figure 16 (bottom panel) 
by summing the column densities of all compo- 
nents per system as given in their Tables 2-10. 
If there were components with upper limits for 
column densities, we set the total system column 
density to an upper limit if the components with 
upper limits were more than 30% of the total. If 
there were components with lower limits for col- 
umn densities, we set the total column density 
to a lower limit. There was no case when these 
criteria conflicted. For the high-redshift sample, 
there were 39 detections and one upper limit for 
N(Si +3 )/N(C+ 3 ), for doublets with the observed 
low-redshift column density limits. 



Since both high- and low-redshift samples con- 
tained at least one upper limit, we used survival 
analysis to enable those limits to contribute statis- 
tically. We used the Astronomy SURVival Anal- 
ysis package (ASURV Rev. 1.3, last described in 



iLavallev et aL1ll992h to compare the two iV(Si +3 )/ 



N(C +3 ) data sets. First, we tested whether the 
low-redshift A r (Si +3 )/A(C+ 3 ) distribution shared 
the same parent population as the high-redshift ra- 
tios. From several univariate ASURV statisticsF^I 
we conclude that the two populations are statis- 
tically similar (i.e., the null hypothesis cannot be 
ruled out with high confidence). 

Next we measured the median ratio of the 
parent population with the Kaplan-Meier estima- 
tor The estimated median of the combined 
low- and high-redshift data sets was (A(Si +3 )/ 
A(C+ 3 )) = 0.16, and the 25th and 75th per- 



12 For more information about univariate analyses used here 
(the two Gehan' s, the Peto-Peto, and Peto -Prentice gener- 
alized tests), see [Fcigclson & Nelson] l ll985f) . 

13 For a useful description of the Kaplan -Meier estima t or in a 
context similar to that used here, see lSimcoe et al. 
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centiles were 0.09 and 0.26, respectively. 

Though the estimated means of the high- and 
low-redshift samples indicated that there should 
be no evolution of the ratio with redshift, we 
checked for a correlation^ Once again, the null 
hypothesis (i.e., that there is no correlation) can- 
not be ruled out with high confidence. 

4.5. Nature of Systems with Si IV and 
C IV Absorption 

We concluded that there has been no evolution 
in V(Si +3 )/V(C+ 3 ) from z = 4.5 -»• 0, based 



on ou r sample and that from iBoksenberg et al 
(2003j). Next we explored what the lack of evo- 



lution means, and we began by disentangling the 
physics involved in the ratio V(Si +3 )/A(C+ 3 ): 



V(Si 4 



N(C+ 3 ) 




(5) 



where Lx is the size of the cloud (enriched with X) , 
along the line of sight; nx is the volume density 
of element X; and Xx i s * ne fraction of X ionized 
into ion X'. The first term on the right-hand side 
is affected by the spatial distribution; the second, 
the metal abundances; and the third, the ionizing 
background. 

The Universal trend is for structure to collapse 
and become denser with age in a ACDM Universe 
and for feedback processes to disperse and mix 
metals on varying scales. Therefore, more of the 
filamentary structure is enriched as the Universe 
ages; though feedback may preferentially enrich 
voids instead o f filaments, as a result of the den- 
sity difference ( Kawata fc Rauch 20071 ). However, 
in our systems with both Si IV and C IV absorp- 
tion, the absorption profiles trace each other quite 
well. There are no obvious system where 7V(Si ) / 
A(C +3 ) varies significantly between the compo- 
nents. Therefore, we infer Lsi ~ Lq- 

We know that the metallicity of the Universe in- 
creases with age, but the relative abundance of sil- 
icon and carbon does not, necessarily, follow suit. 
Hence, the evolution in nsi/nc is unclear. 



14 For more information about the bivariate (Cox propor- 
tional hazard model, generalized Kendall's tau, and Spear- 
man's rho) and linear regression analyses (E M algorithm 
and B uckley-James method) used here, see llsobe et alj 
Jl986h . 



We explored the effect of the ionizing back- 
ground with a suite of simple photoionization 
models, in order to develop our understanding of 
the last term in Equation [5] We used the spec- 
tral synth esis program CLOUD Y v08, as last de- 
scribed by iFerland et all (1998). We modeled the 
medium as a plane-par allel slab, ionized by the 
Haardt fc Madaul (|l996l updated in 2005) ultravi- 



olet background, for quasar (Q), galaxy (G), and 
quasar+galaxy (Q+G) models. We set the num- 
ber density of hydrogen tin = 0.1cm -3 , though 
our models were insensitive to this parameter in 
the optically-thin regime. We assumed a neu- 
tral column density logiVni = 16 and metallicity 
Z = 0.001 Zq. We tested two cases: solar relative 
abundances and Si-enhanced; the increase in sil- 
ic on was such that [Si / C] = +0.77, as measured 
in lAguirre et al.l ((2004). We varied the ionization 
parameter log U = [—5, —4, —3, —2, —1], which is a 
dimensionless ratio of the flux of hydrogen-ionizing 
photons to nn ■ c. 

All models could reproduce the observed lack 
of evolution in the ionic ratios with redshift. In 
Figure El we show V(Si +3 )/V(C+ 3 ) as a func- 
tion of redshift for the three UVB models, with 
solar relative abundances and Si-enhanced, and 
for \ogU = —2. We did not fit the observa- 
tions, but they are reproduced well by the solar 
relative abundance, Q+G model. More impor- 
tantly, all UVB models could reproduce a shape 
consistent with the observed lack of evolution in 
N(Si +3 )/N(C +3 ), given the right choices of e.g., 
logCZ, logAni, [Si/C]. The overall magnitude is 
nearly freely scalable by adjusting these param- 
eters, because the observed sample of Si IV and 
C IV doublets are not drawn from a single type of 
cloud, as we modeled. Evidently, there is no need 
for a particularly soft UVB (i.e., model G) to re- 
produce the lack of redshift evolution in V(Si +3 )/ 
N(C+ 3 ). 

For there to be no evolution in V(Si +3 )/ 
A"(C +3 ) from z — 4.5 — >• 0, the abundance, 
ionizing background, and structure of silicon- 
and carbon-enriched gas are constrained to be in 
"balance." The observations indicate that these 
three processes evolved to maintain a constant 
ratio of (V(Si +3 )/V(C+ 3 )) = 0.16 for nearly 
12Gyr, for absorbers with logA r (Si +3 ) > 11.9 
and logV(C+ 3 ) > 13.37. 

Disentangling the effect of the detailed physics 
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Fig. 8.— Ionic ratio V(Si +3 )/iV(C +3 ) as a function of V(C+ 3 ) (left) and Z ] *<n (right). The 2 < z^ a < 4.5 

Tables 
High- 



data (open circles) come from the total column densities per system from iBoksenberg et al. 
2-10). The low-redshift data have logiV(Si +3 ) > 11.9 and log V(C+ 3 ) > 13.37 (black, filled squares) 
redshift systems meeting these limits are the larger, black circles, while the rest of the sample are the smaller, 
gray circles. Survival analysis shows that the high- and low-redshift ratios are drawn from the same parent 
population and do not evolve significantly with redshift. 



(e.g., changing silicon and carbon abundances, 
variation in the physical properties of absorbing 
clouds) require cosmological hydrodynamic simu- 
lations that could resolve enrichment processes in 
galactic halos and the large-scale structure. 

5. Summary 

We conducted a blind survey for z < 1 Si IV 
doublets in the HST UV spectra of 49 quasars. We 
identified 22 definite Si IV systems (G = 1) and six 
"high- likely" ones (G — 2), and this represents the 
largest sample of low-redshift Si IV doublets prior 
to Servicing Mission 4. From a sample of 20 z < 1 
Si IV doublets with both lines detected at W r > 
3aw r , we measured a line density dA/"siiv/dX = 
1.4±8:a for logiV(Si +3 ) > 12.9. 

We constructed frequency distributions of the 
column densities /(V(Si +3 )) and the rest equiv- 
alent widths /(W rj i393). Both were approxi- 
mated well by power laws. The best-fit power 
law to /(7V(Si+ 3 )) had slope a N = — 1.61±g:|f 
and normalization k = (1.18^03(5) x 10 _14 cm 2 . 
We compared the Si IV line density to high- 
redshift observations by integrating f(N(Si +3 )), 
and cbVsiiv/dX does not evolve significantly from 



From the first moment of /(iV(Si +3 )), we mea- 
sured the Si +3 mass density relative to the criti- 
cal density: n si+3 = (3.71±?;||) x 10~ 8 for 13 < 
logiV(Si+ 3 ) < 15. This value was estimated, with 
Monte Carlo sampling of the distributions, to be 
a factor of 4.8^ 3 'g greater than th e measur e ments 



fro m the 2 < z < 4.5 studies from lSongailal ([2001) 
and IScannapieco et al.1 (|2006h . 

From a simple linear fit, we estimated the rate 
of increase in the Si +3 mass density over time to be 
dft Si+ 3/dt agc = (0.61±0.23)xlO- 8 Gyr. Though a 
linear model is extremely simplistic and not phys- 
ically motivated, it does lend support to the z <1 
figi+3 being a true increase over the 1.5 < z < 5.5 
observations, which, when fit by themselves, result 
in no statistically significant temporal evolution. 

Any increase in ^gi+ 3 i s probably driven by 
the increase in the number of high-column density 
absorbers (i.e., shallower ajv), since dJ\fsnv/dX 
does not increase significantly from z«3->0. 
We also compared the ionic ratio V"(Si +3 )/ 



iV(C +3 ) from the current st udy and IPaper I 



with t he high-redshift sample of lBoksenberg et al 
(l2003h . From survival analysis of the two pop- 
ulations, we concluded that the ionic ratios of 
the high- and low-redshift distributions are drawn 



16 




-1393 



Fig. 9. — Predictions from simple photoioniza- 
tion models for A(Si +3 )/A(C+ 3 ) as a function of 
Z1393. All models could reproduce the observed 
lack of evolution in V(Si +3 )/A(C+ 3 ) with redshift 
and the overall mag nitude of N(Si +3 )/N(C+ 3 ) 
through careful choice of the CLOUDY parame- 
ters (e.g., log £7, [Si/C]). For this plot, the ion- 
ization parameter is \ogU = —2; the metallicity 
is 0.001 Zq. The solid lines indicate models with 
solar relative abundances and the dashed lines, Si- 
enhanced (i.e., [Si/C] = +0.77). The gray, back- 
ground points are the observations (black symbols 
from Figure [SJ right), but the CLOUDY models 
were not fit to the observations. 



from the same parent population, with median 
(A(Si +3 )/A(C+ 3 )) = 0.16. The lack of evolution 
in A(Si +3 )/A(C+ 3 ) from z = 4.5 -> places con- 
straints on the evolution of the metal production, 
feedback processes, and the ionizing background. 
These three processes evolve in some balanced 
fashion to have Si IV and C IV absorbers evolve 
in lock-step for the last w 12 Gyr. 

We explored the effect of the ionizing back- 
ground on the (non)evolution of A(Si +3 )/A(C+ 3 ) 
with a suite of simple CLOUDY models. We 
varied the backg r ound by using the canonical 
Haardt fc Madau (1996) quasar, galaxy, and 
quasar+galaxy models. All three backgrounds 
could result in relatively constant N(Si )/ 



N(C +3 ) over redshift, for models with the "right" 
choices of e.g., ionizing parameter, [Si/C]. There- 
fore, a soft UVB (i.e., model G) is not preferred. 
In general, more observations — at low and high 



redshift — are needed to increase the statistical sig- 
nificance of the trends that are currently highly 
suggestive. We eagerly anticipate new low-redshift 
results from the HST Cosmic Origins Spectro- 
graph fCOS. iMorse et al.|[l998h . Meanwhile, cos- 
mological hydrodynamic simulations should be 
leveraged to understand how metal production 
and dispersal and the ionizing background inter- 
act to evolve Si IV and C IV absorbers in tandem. 
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